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We study the squeezing as well as the statistical properties of the output radiation from a de- 
generate parametric oscillator coupled to a squeezed vacuum reservoir employing the stochastic 
differential equations associated with the normal ordering. It is found that the degree of squeezing 
of the output radiation is less than the corresponding cavity radiation. However, for output radia- 
tion the correlation of the quadrature operators evaluated at different times also exhibits squeezing, 
which is the reason for quenching of the overall noise in one of the quadrature components of the 
squeezing spectrum even when the oscillator is coupled to a vacuum reservoir. Moreover, coupling 
the oscillator to the squeezed vacuum reservoir enhances the squeezing exponentially and it also 
increases the mean photon number. 

PACS numbers: 42.50.Dv, 42.65.Yj, 42.50.Ar 



I. INTRODUCTION 

Optical degenerate parametric oscillator is one of the 
most interesting and well studied devices in the nonlin- 
ear quantum optics p], [H, H, 0, IE @> • I n a degenerate 
parametric oscillator, a pump photon of frequency 2u> is 
down converted by nonlinear crystal into a pair of sig- 
nal photons each of frequency to. Due to this inherent 
two-photon nature of the interaction, the parametric os- 
cillator is found to be a good source of a squeezed light. 
The maximum degree of squeezing of the cavity radiation 
coupled to an ordinary vacuum is found to be 50% at the 
critical point by many authors using various approaches 
0, 0- S] • This limitation in squeezing is associated with 
the leakage through the mirror and the inevitable ampli- 
fication of the quantum fluctuations in the cavity due to 
thermal heating. Since the vacuum field has no definite 
phase, squeezing of the cavity radiation when coupled 
to a vacuum reservoir would be degraded. However, if 
the ordinary vacuum is replaced by a squeezed vacuum 
reservoir, it is expected that the fluctuations entering the 
cavity is biased and hence the squeezing of the cavity ra- 
diation could be enhanced @, U 0, H[ provided that the 
reservoir is squeezed in the right quadrature. 

Theoretical analysis of the quantum fluctuations and 
photon statistics of the cavity radiation of a degener- 
ate parametric oscillator coupled to a squeezed vacuum 
reservoir has been exhaustively made by many authors 
P) B El ) although it is the output radiation which is prac- 
tically accessible for an experiment. To our knowledge, 
there is no a thorough study of the squeezing and statis- 
tical properties of the output radiation except for some 
discussions related with the squeezing spectrum [H, H, @ . 
For instance, Collett and Gardiner [f| claimed, after com- 
paring the variance of the quadrature operators for the 
cavity radiation with the squeezing spectrum for the out- 
put radiation, that the squeezing of the output radiation 



is much better than the cavity radiation. Since the out- 
put squeezing spectrum corresponds to the correlation of 
the output quadrature operators evaluated at different 
times, t and t + r, in the frequency domain, we believe 
that it would not be right comparing the squeezing spec- 
trum with the variance of the cavity radiation. 

We hence devoted this work to the analysis of the 
output radiation of a degenerate parametric oscillator 
coupled to a squeezed vacuum reservoir using the input- 
output relation introduced by Gardiner and Collett 
We device a means of evaluating a correlation between 
the noise forces arising from the reservoir and vacuum 
fluctuations in the cavity that required in studying the 
properties of the radiation outside the cavity. Moreover, 
we followed an adhoc approach in casting the charac- 
teristic function describing the cavity radiation onto the 
corresponding output variables. On the basis of these ap- 
proaches we study the squeezing as well as the statistical 
properties of the output radiation employing the stochas- 
tic differential equations associated with the normal or- 
dering. In particular, we calculate the quadrature vari- 
ance, squeezing spectrum, photon number distribution, 
and power spectrum for the output radiation. We also 
compare the results we obtained with the corresponding 
values for the cavity radiation. We believe that the ap- 
proach followed in this work could also be used to study 
the output radiation of other quantum optical systems. 

II. STOCHASTIC DIFFERENTIAL 
EQUATIONS 

The degenerate parametric oscillator can be described 
in the interaction picture and in the rotating-wave ap- 
proximation, upon treating the pump radiation classi- 
cally, by the Hamiltonian of the form 

£ = f(« t2 -& 2 ), (1) 
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where e is proportional to the amplitude of the external 
coherent radiation which is taken to be a real-positive 
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constant and a is the annihilation operator for the signal 
mode. We consider the case in which a continuum mode 
of squeezed vacuum centered at frequency w is allowed to 
enter the cavity through one of the coupler mirrors. In 
this case, the master equation describing the degenerate 
parametric oscillator coupled to a broadband squeezed 
vacuum reservoir in the interaction picture is found using 
the standard procedure [10( to be 



- = -^ p-ap-pJ +pa 
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where N and M represent the squeezed vacuum reservoir, 
with N = sinh 2 (r), M — sinh(r) cosh(r), k is the cavity 
damping constant, and r is the squeeze parameter. 

We next proceed to obtain the stochastic differential 
equations pertinent to the cavity mode variables using 
the master equation @ . To this end, employing the fact 
that 



it is possible see that 



-(a(t))=--(a(t))+e(aHt)), 



-(a\t)a(t)) = - K (a\t)a(t)) 



(3) 



(4) 



+ e[(a 2 (t)) +(a t (t))]+KN, (5) 



— (a 2 (t)) = -n(a 2 {t)) + 2e(aUt)a(t)) +e + nM. (6) 
at 

We notice that the operators in Eqs. |@J, ||SJ|, and ([6]) 
are in the normal order. Hence the corresponding ex- 
pressions in the c-number variables associated with the 
normal ordering take the form 



-~(a{t))+e{a*{t)), 



(7) 



-(a*(t)a(t)) = - K {a*(t)a(t)) 



+ e[(a?(t)) + (a* (t))]+KN, (8) 



— (a 2 (t)) = -n(a 2 (t)) + 2e(a* (t)a(t)) +e + kM. (9) 
at 



On the basis of Eq. (JT]) , one can write that 

j t a(t) = -^a(t) + ea*(t) + F(t), (10) 

where F(t) is a complex Gaussian white noise force the 
mean and correlation properties of which remain to be 
determined. It is not difficult to see that Eq. ([7|) will be 
equal to the expectation value of Eq. (TU)) provided that 

(F(t)) = 0. (11) 

On the other hand, one can verify applying Eqs. (0, {(5J), 
©, and (|10p along with the fact that the noise force at 
time t does not correlate with the cavity mode variables 
at earlier times that 

(F(t')F(t)) = (KM + e)S(t-t'), (12) 

(F*(t)F{t)) = KNS(t - t'). (13) 

We observe that Eqs. ([IT]). l[T2]). and (T3]) represent the 
mean and correlation properties of the noise force. 

In general, it is possible to express Eq. (fl7J|) in a more 
convenient from, upon introducing two variables defined 
by 



as 



with 



a±(t) = a*(t) ± a(t), 



^a± = -)^\±a± + F±(i), 



(14) 
(15) 



A±=k=f2£, (16) 

F ± (t) = F*(t)±F(t), (17) 
where the formal solution of Eq. (|15p can be written as 

a±(t +r) = a±(t)e 2 T 

e-^ {T - T,) F±{t + T ')dT' . (18) 



Now in view of Eqs. fM]) and (| 18|) into account, we get 
a(t + r) = A + (r)a(t) + A_(r)a*{t) 

+ B + (t + T)~ fl_(t + T), (19) 

in which 



A±(r) = i[ e -^±e-^], 



(20) 



B±(t + r) = i [ e-^ T - T ^F ± {t + T')dT'. (21) 
2 Jo 

Next upon setting i = in Eq. (fT§|) and then replacing 
r by i, we have 

a(i) = A+(t)a(0) + A-(t)a*(0) + B+(t) - B-(t). (22) 

It perhaps worth mentioning that Eqs. (fT9")> and (|2"2")) is 
applied in calculating various quantities of interest. We 
also realize that these solutions would be well-behaved 
functions at steady state, if k > 2s. Hence we designate 
k = 2s as a critical point. 
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III. QUADRATURE FLUCTUATIONS 

The squeezing properties of a single-mode output ra- 
diation can be described with the aid of the quadrature 
operators defined by 



+ a a 



and 



= K&out - O-out). 



(23) 



(24) 



A. Quadrature variance 

The variances of the quadrature operators (|23p and 
(|24|) can be put, using the boson commutation relations 
for the output radiation, in the form 

Aa 2 ±(out) = 1 ± [<&£*(*)) + ± 2(aLt(t)a ut(t)) 

+ (&Lt(t)f + {a ut(t)) 2 ± 2{al ut (t)){a out (t))} . 

(25) 

We notice that the operators in Eq. (j25|) are in the nor- 
mal order. Hence the corresponding expression in terms 
of the c-number variables associated with the normal or- 
dering would be 

Aa 2 ±(out) = 1 ± [<a£ t (i)) + (a 2 out (t)) ± 2(a* out (t)a out (t)) 

+ «uM) 2 + (<w(i)> 2 ± 2{a* out {t)){a out {t))} . 

(26) 

It is a well established fact that the output radiation can 
be defined in terms of the cavity mode variables 0] as 



t (t) = vM*) - -p-pflW. 



(27) 



where F R (t) is the noise force associated with the reser- 
voir modes and satisfies the correlations: 



(Mm = o, 

(F R (t)F R (t')) = KN6(t-t% 
(F R (t)F R (t)) = K M6(t-t'). 



(28) 
(29) 



(30) 



If the cavity mode is taken to be initially in a vacuum 
state, Eq. (|2"6")l reduces to 



Aa 



= l + 2{a* out (t)a out {t)) 
± [« 2 ut(t)) + (a 2 out (t))]. 



(31) 



We next seek to determine the correlations involved in 
Eq. (f3"Tj) . To this end, making use of Eq. l[2"7]). we see 
that 

(a* out (t)a out (t)) = K(a*(t)a(t)) + ~(F R (t)F R (t)) 

-{a*{t)F R (t))-{F R (t)a(t)). (32) 



We observe that the noise force F(t) represents the con- 
tribution from the cavity vacuum fluctuations and reser- 
voir, and hence can be expressed 



F(t) = F c (t) + F R (t), 



(33) 



where Fc(t) is the noise force corresponding to the sys- 
tem in the cavity in the absence of damping. Since the 
noise force for the reservoir F R (t) does not correlate with 
Fc(t) and the system variables at the earlier times, it is 
not difficult to verify that 

(F R (t)a*(t)) + (F R (t)a(t)) = kN. (34) 

Thus on account of Eqs. (J29J), (J32J) , and (J34J) , we get 

(a* out (t)a out (t)) = K(a*(t)a(t)) + N(l - «). (35) 

We are usually interested in steady state values in quan- 
tum optics. In this regard, it is easy to check applying 
Eqs. © and © that 



(a*(t)a(t)) ss = 



2e 2 



k(Nk + 2eM) 



(36) 



as a result the mean photon number for the output radi- 
ation at steady state takes the form 



t(t)a ou t{t)) t 



N + 



2ke 2 



2en(2Ne + kM) 



Ae 2 



k? - 4e 2 



(37) 
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FIG. 1: Plots of the mean photon number of the output and 
cavity radiation at steady state for k = 0.8 and r — 0.75. 



We notice that the first term in Eq. (f3"T)) is the contri- 
bution from the squeezed vacuum radiation that exists 
outside the cavity, the second term corresponds to the 
light produced by the parametric oscillator and then es- 
capes through the mirror, whereas the last term repre- 
sents the correlation between the oscillator and reservoir 
due to the associated noises. We realize that the larger 
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r=0.75 
M).5 
r=0.25 





FIG. 2: Plots of the mean photon number of the output ra- FIG. 3: Plots of the quadrature variance for the output and 
diation at steady state for k = 0.8 and different values of cavity radiation at steady state for n — 0.8 and r — 0.75. 



the intensity of the external coherent radiation, the more 
probable it is converted into the cavity radiation by the 
crystal. It is a known fact that the squeezed parame- 
ter is directly related to the mean photon number of the 
squeezed vacuum modes. Thus the more the reservoir 
modes are reflected from the mirror, the more the mean 
photon number outside the cavity would be. We hence 
clearly see from Fig. 2 that the mean photon number for 
the output radiation increases with the amplitude of the 
coherent radiation and squeeze parameter of the reser- 
voir. However, in spite of the presence of the squeezed 
vacuum radiation outside the cavity, we clearly see from 
Fig. 1 that the mean photon number is greater in the 
cavity, since only few of the photons in the cavity can 
escape through the mirror. 

Furthermore, it is possible to verify following a similar 
procedure that 



(a 2 out (t)) ss =M + 



k 2 e 2ke(kN + 2eM) 



k - 4e 2 



k - 4e 2 



(38) 



so that the variance of the output quadrature operators 
are found to be 



A 4(ou f ) 



,±2r 



1± 



2en 



k T 2e 



which reduces at critical point to 



Aa +(out) = °°> 



Aa 2 _ (out) = e~ 2r [1 - 0.5fc] 



(39) 



(40) 



(41) 



In addition, the quadrature variances for the cavity radi- 
ation turn out following a similar approach to be 



Aal = , 



,±2r 



(42) 



in which at critical point Aa^_ — > oo and Aa_ = 



l D ±2r 



Since k < 1, one can easily infer from Eq. (|41j) that the 
squeezing of the output radiation is less than 50% for a 
vacuum reservoir at critical point. We also see from Fig. 
(3) that the degree of squeezing of the cavity radiation is 
far better than the output radiation except for very small 
values of the amplitude of the driving radiation. The 
unbiased noise fluctuations arising while the radiation 
crosses the boundary (coupler mirror) is believed to be 
the reason for the reduction of the degree of squeezing 
for the output radiation. 



B. Squeezing spectrum 

In this section, we seek to calculate the squeezing spec- 
trum for the output radiation, which can be defined as 



S£ t (w) = 2Re e wr {a°^{t + T), a°^{t)) ss dT, (43) 
Jo 

where 

-(a%*(t + T)){dg*(t)). (44) 

With the application of (O, {23]), and Eq. g3D 
can be put, for the cavity radiation initially in vacuum 
state, in terms of the c-number variables associated with 
the normal ordering as 

poo 

S^ ut (co) = l±2Re / e^ia^it + T)a ^\t)) ss dT. 
Jo 

(45) 

Now making use of Eq. (|27|) , we write 

+ T)a°^(t)) ss = K(a±(t + r)a±(t))„ 

+ -(F R ±(t + T)F R± {t)) ss 

K 

~ (a±(t + T)F R ±(t)) ss 
-{a±(t)F R ±{t + T)) ss , (46) 
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in which 



IV. PHOTON STATISTICS 



F R± (t)=F^(t)±F R (t). (47) 

On the basis of the fact that the noise force does not cor- 
relate with system variables at earlier times, it is possible 
to verify that 



(a±(t + T)a±(t)) ss = 



k T 2e 



[k(M±N) 



■ E\e 2 



(48) 



(F R ±(t)F R± (t + T)} = 2k(M±N)S(t), (49) 

(a±(t + T)F R ±(t)) =2 K (M±N)e~ t r T , (50) 

as a result employing Eqs. (|15j). (14^) . (gHJ) , (gH]), and 
([SUf. we find 

8£K 



„±2r 



1± 



(k t 2e) 2 + 4w 2 
At critical point and w = 0, we see that 5°"'(0) = 




FIG. 4: Plots of the squeezing spectrum of the output radia- 
tion at steady state for k — 0.8, r = 0.75, and different values 
of e. 



We notice from Eq. (|?T|) that for e = 0, the squeezing 
spectrum is completely flat with its value depending on 
the squeeze parameter. However, for s 7^ the depth 
of the squeezing spectrum at to = increases with the 
amplitude of the coherent radiation and it would be zero 
at critical point irrespective of the values of the squeeze 
parameter. It is not difficult to see that the width of the 
squeezing spectrum deceases with the amplitude of the 
coherent radiation. We realize upon comparing the re- 
sults shown in Figs. 3 and 4 that the squeezing spectrum 
can take smaller values than the corresponding quadra- 
ture variance of the output radiation for the same r and 
e. But this does not mean that the squeezing of the out- 
put radiation is greater than the cavity radiation, since 
the squeezing spectrum corresponds to the correlation of 
the quadrature operators evaluated at different times un- 
like the quadrature variance of the cavity radiation which 
is associated with the correlation of the quadrature vari- 
ance at equal time. 



A. Photon number distribution 

The probability for finding n photons in a single-mode 
radiation can be expressed in terms of the pertinent Q- 
function as 



P(n) 



7T d 



2n 



n\ da n da* 



-[Q(a)exp(a*a)] Q , =Q=0 , (52) 



where the Q-function for the single-mode radiation can 
be defined as 



Q(M) 



1 



d 2 z 4>{z, z* , t) exp[z*a — za*], (53) 



in which the antinormally ordered characteristic function 
(51) (j)(z,z*,t) is given by 

<j)(z, z*, t) = Tr{p(0)e- z * a(t) e zat W}. (54) 
Using the identity 

eV = e¥e'^l (55) 

Eq. (|54p can be written in terms of the c-number vari- 
ables associated with the normal ordering as 



(z,z*,t) = e- z ' z (exp(za*(t) - z*a(t))). 



(56) 



We note that a is a Gaussian variable with zero mean. 
Hence one can readily verify that [TTI ] 



(exp(za*(i) — z*a(i))) = exp 



\{{za*{t)-z*a{t)f) 



(57) 



Therefore, it is possible to express the antinormally or- 
dered characteristic function for the output radiation as 



e 2 z exp 



z 2 <cC t (t))+z* (a 2 out (t)) 



-2z*z{a* out (t)aout(t)))}, (58) 
which can also be put at steady state in the form 

<j) out (z, z*) = exp[-z*za + ^(z 2 + z* 2 )}, (59) 

where a = 1 + h and b = (a 2 ut (t)) ss , in which n repre- 
sents the mean photon number for the output radiation 
at steady state (|3"7) . Hence substituting Eq. ([55]) into 
([53")) and then performing the integration, the Q-function 
describing the output radiation is found to have the form 



S («) 



Vu 2 - « 2 



■ exp 



-ua a + — (a + a ) 



(60) 



with u = 3 a .-i and v = 

a 2 —b 2 



b 

a 2 -b 2 ■ 
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On the other hand, insertion of Eqs. (|60)) into ([52 
results in 



P out (n) 



2 



cxp 



da* n da n 
(1 - u)a*a + ^(«* 2 + a 2 ) 



a—a*—0 

(61) 



from which follows 



P out (n) 



X 



(1 - u) l v k+3 



2 k +H\j\k\ 

i,j,k=0 



{2k + i)\ (2j + z)! 
(2k+i-n)\ (2j+i-n)\ 



xa 



2j+i-n 



a— a* —0 



(62) 



Now applying the condition, a* = a = 0, we note that 
Eq. (I62|) would be different from zero provided that 2k + 
i — n = and 2j + i — n = 0, as a result k = j and 
j = ^y 1 . Consequently, the photon number distribution 
for the output radiation finally takes the form 



P out {n) = 



nl 



[l + 2n + n 2 -(a 2 out m 
[n 2 +n~ (a 2 out (t)) 2 ss . 



t (*)>?r 



i=0 



2»-*il[(^i)!] 2 



(63) 



We notice that there is a finite probability for counting 
odd number of photons outside the cavity, even though 
degenerate parametric oscillator generates pairs of pho- 
tons. This could be related to the fact that odd number 
of photons can escape through the coupler mirror. 



B. Power spectrum 



as a result the power spectrum (164[) for the output radi- 
ation turns out to be 



S ou \lu) = N + 2ke 



1-2N+ 2M 



1 + 2N + 2M 

(k - 2s) 2 + 4w 2 + (k + 2e) 2 + 4oj 2 



(66) 



In a similar manner the power spectrum for the cavity 
radiation is found to be 



k(N + M) + e k(N -M)-e 



(k - 2e) 2 + Alo 2 (k + 2e) 2 + Au 2 



(67) 




FIG. 5: Plots of the power spectrum of the output and cavity 
radiation at steady state for k = 0.8, r — 0.5, and e = 0.2. 



We realize from Eqs. (|66[) and (j67|) that the width of the 
cavity and output power spectra is independent of the 
squeeze parameter, but it decreases with the amplitude 
of the driving radiation. As clearly shown in Fig. 5 the 
power spectrum for cavity radiation is more picked and 
narrowed. However, the width for the cavity and output 
radiation is the same for the same values of e. 



The power spectrum for the output radiation can be 
expressed in terms of the c-number variables associated 
with the normal ordering as 

/>oo 

S OM *M = 2i?e / (a* out {t)a out (t + r)) ss e MT . (64) 
Jo 

Here with the aid of Eqs. Q2J), {T3}, (JTSJ), (2Q), (HQ), 
([27]). ([%?])■ and §U§, we see that 

(a* out {t)a out (t + T)) ss = |((a*(t)a(t)) ss + (a 2 (i)) ss 

- (M + N))e-^ T 

+ ^({a*(t)a(t)) ss -{a 2 (t)) ss 

— (N — M))e-^ T + NS(t), 

(65) 



V. CONCLUSION 

In this paper, we present a thorough study of the 
squeezing and statistical properties of the output radi- 
ation generated by a degenerate parametric oscillator 
coupled to a squeezed vacuum reservoir using the input- 
output relation introduced by Gardiner and Collett 
We also make a comparison with the results of the cavity 
radiation. For both the cavity and output radiations a 
maximum squeezing is occured at a critical point, k = 2e. 
The maximum squeezing for a cavity radiation when the 
degenerate parametric oscillator is coupled to a vacuum 
reservoir is found to be 50%, irrespective of the cavity 
damping constant. The same result has been obtained 
by many authors [1-4] using various approaches. How- 
ever, we see from Eq. (|4ip that the squeezing for the 
output radiation increases with the damping constant at 
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critical point. One can also clearly see from Fig. 3 that 
the squeezing of the cavity radiation is greater than that 
of the output radiation. 

Coupling the degenerate parametric oscillator with the 
squeezed vacuum reservoir is found to exponentially en- 
hance the squeezing of the output as well as the cavity 
radiation. Our calculation of the squeezing spectrum for 
the output radiation indicates that a complete quench- 
ing of the noise in one of the quadrature components is 
possible at the critical point, even when the degenerate 
parametric oscillator is coupled to the vacuum reservoir. 
Since the squeezing spectrum corresponds to the corre- 
lation of the quadrature operators at different times, we 
note that the complete suppression of the noise could be 
achieved if correlation at different times results squeezing 
and there is a way of measuring the overall squeezing over 
sufficiently long period of time. We also see from Fig. 4 
that the width of the squeezing spectrum decreases with 



the amplitude of the driving radiation. 

As the degree of squeezing of the output and cavity 
radiations increases with the amplitude of the coherent 
radiation and squeeze parameter so does the mean pho- 
ton number. This indicates that the degenerate para- 
metric oscillator produces not only a light with high de- 
gree of squeezing but also a quite intense light at the 
critical point. Moreover, even though only even number 
of photons are produced by the degenerate parametric 
oscillator, there is a finite probability for counting odd 
number of photons outside the cavity. Furthermore, the 
power spectrum for the cavity and output radiations is 
found to be picked at uj — 0. As clearly shown in Fig. 
5 the height of the power spectrum for the cavity radia- 
tion is greater than that of the output radiation, and yet 
have the same width. It is not difficult to see from Eqs. 
(l66|) and (|6T[) that the width of the power spectrum is 
independent of r. 
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